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INTRODUCTION 


Classical  energy  and  dynamical  budgeting  studies  have 
traditionally  encompassed  global  or  near  global  regions. 

Many  more  recent  studies,  however,  have  shown  that  limited 
area  budgets  can  provide  some  excellent  insight  into  phenom¬ 
ena  of  limited  geographical  size  if: 

a)  the  boundaries  of  the  system  can  be  defined,  and 

b)  the  horizontal  flux  through  the  boundaries  can  be  deter- 
mined . 

While  the  tropics  have  provided  the  greatest  poten¬ 
tial  for  limited  area  budgeting,  in  this  and  past  work 
(McGinley ,  1973),  we  consider  systems  in  the  westerlies  where 
the  boundary  flux  problem  can  be  significant.  This  problem 
is  diminished  somewhat  by  the  use  of  high  resolution  varia¬ 
tional  analysis  schemes. 

This  present  work  is  a  continuation  and  offshoot  of 
developmental  work  on  a  three  dimensional  analysis  procedure 
for  mountainous  terrain,  including  analysis  of  basic  meteoro¬ 
logical  parameters  and  computation  of  terrain  induced  and  mass 
constrained  vertical  velocity  (Sasaki,  et  al.,  1977).  The 

fine  tuning  of  this  scheme  continued  into  the  present  contract 
period  . 


Modifications  and  improvements  to  the  basic  analysis 
included : 

a)  Addition  of  a  strong  inequality  constraint  (see  Appendix) 
to  adjust  super  adiabatic  layers  arising  from  the  analy¬ 
sis  (a  manuscript  is  being  prepared  for  publication) . 

b)  A  routine  to  extract  vertical  cross  sections  from  the 
sigma  system.  This  facilitated  case  study  analysis. 

Modifications  not  completed  but  under  way  are: 

a)  Merging  a  first  guess  field  into  the  analysis  scheme  for 
operational  use  or  for  extended  case  studies. 

b)  A  scheme  to  remove  obviously  bad  radiosonde  data  prior 
to  objective  analysis. 

The  budget  analysis  considered  bulk  or  volume  evalu¬ 
ations  of  fundamental  energy  and  dynamical  quantities, 
including  the  various  forms  of  energy,  vorticity,  potential 
vorticity  and  frontal  intensity.  Work  by  Ninomiya  and 
Matsumoto  (1968,  1969)  was  used  as  a  basis  for  this  work. 

In  addition,  energy  and  dynamical  tendencies  evaluated 
from  equations  for  kinetic,  internal,  latent  heat,  vorticity 
were  computed  and  compared  with  the  bulk  energetics  and 
dynamics . 

The  case  studies  include: 

a)  6  November  1973,  strong  cold  front  crossing  Alps  into 
Mediterranean  Sea . 

b)  21/22  July  1976,  strong  summer  cold  front,  over  Alps  into 
Mediterranean . 
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This  and  succeeding  work  should  serve  as  a  benchmark 
for  numerical  studies  conducted  by  others  (Buzzi  and  Tribaldi, 
1977;  Trevisan,  1976;  and  Bleck,  1977).  The  problem  of 
Mediterranean  cyclogenesis  is  a  complex  one  involving  a 
mountain  range  of  sub-synoptic  scale,  a  persistent  thermal 
and  latent  heat  energy  gradient  and  abundant  upper  level 
disturbances  moving  into  the  area  from  the  northwest  to  west. 
The  energetic  and  dynamic  influences  of  these  stationary 
terrain  features  figure  prominently  in  the  consequent  devel¬ 
opment  and  modification  of  moving  weather  systems.  Hopefully 
some  measure  of  these  effects  and  speculations  regarding 
their  origin  can  be  made  within  this  report. 
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BASIC  ANALYSIS 


2  . 

[For  details,  see  (Sasaki,  et  al  .  ,  1977)] 

The  analysis  scheme  uses  twice  a  day  radiosonde 
observations  subjecting  the  data  to  the  following: 

a)  Variational  low  pass  filtering  to  smooth  the  data. 

b)  Adjustment  of  vertical  temperature  profiles  using 
inequality  constraint  (see  Appendix)  . 

c)  Geostrophic  and  hydrostatic  adjustment  of  the  geopoten¬ 
tial  field. 

d)  Calculation  of  surface  pressure  at  the  ground  surface. 

e)  Coordination  of  transformation  from  TT  (=  ln  P  /P)  to  c 

o 

(sigma)  in  the  vertical. 

f)  Application  of  mass  constraint  on  the  wind  field  in  the 
[x,  y,  a]  domain. 

g)  Computation  of  to  (omega)  vertical  motion  fields. 

h)  Computation  of  vorticity  and  potential  vorticity.  ' 

i)  Vertical  cross-sections. 
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3.  BUDGET  CALCULATIONS 


3.1  VOLUME  AVERAGES 

Quantities  are  averaged  over  volume  elements  within 
the  main  three  dimensional  analyzed  domain.  Each  meteoro- 
logical  variable  (A)  is  averaged  in  the  horizontal  and  ver 
tical  according  to: 
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(where  s  is  the  horizontal  area  of  the  volume  at  level  i) 


and 


T 


(where  A  is  the  volume  average  and  O is  the  topmost  O  level). 

Note  that  the  averages  in  the  sigma  system  are  implic¬ 
itly  mass  weighted  since  each  layer  of  fixed  depth  Aa  has  the 


same  mass/per  unit  area,  when  we  consider  some  mean  surface 


pressure  for  the  volume. 

Averages  are  computed  with  n  volumes  and  sub-elements 
for  quantities  of  particular  meteorological  interest,  such  as 
temperature  (enthalpy),  geo-potential  (potential  energy),  wind 
(kinetic  energy),  moisture-mixing  ratio  (latent  heat),  vorti- 
city,  potential  vorticity  and  frontal  intensity  (defined  as 

,  where  n  is  the  horizontal  normal  to  the  frontal  surface)  . 
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3.2  ENERGY  AND  VORTICITY  EQUATIONS 


Equations  for  the  energy  calculations  were  taken  from 
Haltiner  (1971)  and  McGinley  (1973)  and  modified  to  include 
Eularian  and  Lagrangian  forms  in  the  O  system. 

1 •  For  Kinetic  Energy 

a)  Flux  form  (can  be  integrated  to  compute  vertical  flux 
and  can  evaluate  eddy  or  subgrid  fluxes) . 

5T  +  va  •  (KV  +  V  v>  +  h  <K°  +  +  TF  -t?  -  U  -  0 

(where  K  is  kinetic  energy  per  unit  mass,  P  is  surface 

s 

pressure,  $  is  geopotential  of  constant  0  surface,  V  is 

velocity,  w  is  '^)a>  cr  is  The  relation  between  w  (omega) 

and  0  is  (w)a  =  Pg  0  +..  P.  0  .  !Ln  addition,  we  compute  the 

local  tendency  of  kinetic  energy,  |£;  the  Lagrangian  form  — ; 

d  t  d  t 

and  the  parcel  of  budget  volume  form,  where  the  change  is 

defined  with  respect  to  an  arbitrary  moving  element. 

2.  For  Thermal  Energy 

a)  The  Flux  form: 

h  (cppsT>  +  v  •  <VSTV)  +  h  (cpV“>  -  ^ 

where  T  is  temperature,  S  is  a  non  adiabatic  heat  source, 

Cp  is  the  specific  heat  at  constant  pressure,  and  p  is  density. 

We  also  compute: 

9  (C  P  T)  d  (C  P  T)  D(C  P  T) 
p  s  p  s  p  s 

9  t  ’  dt  DT 
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a)  Flux  form: 


3. 


For  Latent  Heat  Energy 


3  (L  qP  s  ) 
3t 


+  V 


<LPs^)  +  3^  <LV°> 


-P  LS 
s 

e 


(where  q  is  mixing  ratio  and  LS  is  a  latent  heat  sink  term) . 

Again,  as  in  1.  and  2.  we  compute 

9  L  P  d  L  P  DLP 
qs  q  s  q  s 

9  t  ’  dt  ’  Dt 


4.  Vertical  Fluxes 

We  can  integrate  Equation  1  to  find  (PsKa)^,  the 
vertical  kinetic  energy  flux.  By  combining  Equations  2  and 
3  and  assuming  that  S  represents  moisture  processes  only, 
we  find  the  vertical  heat  flux: 


OP 


s 


+  C  T ) 

P  CF 


5.  Vorticity  Equations 


Vorticity  (a)  and  potential  vorticity  (b)  are  computed 


as  : 
a)  C 


3v  _  9_u  u  TAN  j>  ,9_v  9_a.  _  9_u  9a. 

3  y  3x  Re  l9a  9  y}  p  9a  9x}PJ 


where  Re  is  earth  radius,  <j>  is  latitude,  and  the  bracketed 
term  is  the  correction  due  to  the  o  system. 


b)  P  =  C 


9_6 
3  P 


_C_  <L0 

ps  3a’ 


0  is  potential  temperature. 


The  vorticity  equation  in  the  pressure  system  is  trans 
formed  into  the  sigma  coordinate  system,  so  that  we  may  eva] 
uate  vorticity  tendency  and  quantify  the  role  of  each  term 
in  vorticity  production. 
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Advection  term 


H  +  »* 

V  (£  +  f)  -  —  -^r - ^  V  •  VP  +  —  Advection  term. 

o  p  ,  3a  s  Pda 

s  s 

+  T~ 

S 

+  (f  +  C)  *  V  -  “  VPg  •  Divergence  term 

s 

,3d)  3u  dcodv  dudcodo*  3v  3(i)  3a^  >.  ...  . 

(s7  ^  -  37  so  +  ^  57JP  •  rs  rs  rjV  Tiltins  term- 
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CASE  STUDIES 


4  . 

As  stated  in  the  Introduction  the  motivation  for  this 
work  is  to  quantify  some  of  the  changes  associated  with  air¬ 
flow  interacting  with  a  mountain  barrier,  specifically  modifi¬ 
cation  to  frontal  boundaries  as  they  pass  over  a  significant 
range . 

Figure  1  shows  the  region  of  analysis  and  data  density 
along  with  terrain  height. 

In  each  time  period  we  will  be  evaluating  changes 
averaged  within  large  volumes  associated  with  surface  and 
upper  level  frontal  zones.  These  volumes  are  positioned 
subjectively  and  are  large  to  insure  reasonable  bulk  evalu¬ 
ations  of  energy  and  dynamical  quantities.  A  similar  tech¬ 
nique  was  used  by  McGinley  (1973).  Figure  2  shows  schemat¬ 
ically  how  the  volume  is  placed  over  a  frontal  surface  moving 
with  the  mean  frontal  velocity  over  the  time  period. 

The  geography  of  the  Mediterranean  region  is  unique 
and  configured  in  such  a  way  that  frontal  zones  either  pass 
over  lpw  hills  in  moving  southeastward  from  the  Atlantic  to 
the  Mediterranean,  or  very  sharp  mountain  ranges  (the  Alps  or 
Pyrenees)  along  adjacent  trajectories.  We  shall  make  use  of 
cross-section  analysis  to  compare  the  "mountain"  versus  "low 
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500  M 

1000-1500  M 
1500-2000  M 
2000-3000  M 


terrain 


DATA  STATION 


Fig.  l. 


Region  of  Analysis,  Rawinsonde  Stations  and  Terrain 
Heights  (m) .  See  key. 
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UPPER  LEVEL(~300mB)  FRONT 


<r=o.i 

~500mB  FRONT 


SURFACE  FRONT 


Placement  of  budget  volume  relative  to  frontal  boundary 
(surface  and  500  mb  fronts  are  illustrated) .  Volume 
propagates  with  movement  of  frontal  surface. 
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terrain"  fronts,  in  addition  to  comparing  their  bulk  proper¬ 
ties.  In  both  case  studies  volume  A  will  be  associated  with 
the  Alps,  volume  B  associated  with  the  Pyrenees,  and  volume 
C  the  terrain  of  southern  France. 

Case  1,  6  November  1973 

Figures  3  and  4  show  frontal  positions  at  surface  and 
middle  levels,  jetstream  positions  and  vertical  motion  fields 
for  00  and  12Z.  The  frontal  surface  moves  over  the  mountain 
ranges  during  this  time  period.  The  frontal  zone  is  divided 
into  three  volume  elements  A  and  B,  which  pass  over  the  Alps 
and  the  Pyrenees  respectively,  and  C  which  moves  over  the 
lower  terrain  between.  Figures  5,  6  and  7  show  changes  of 

energy  and  dynamical  quantities  with  respect  to  these  volumes. 

Note  how  kinetic  energy  in  the  lower  half  of  volumes 
A  and  B  decrease  while  kinetic  energy  in  the  upper  half 
increases.  In  volume  C  we  note  less  kinetic  energy  change 
proportionately. 

Mean  vertical  motion  increases  in  A  and  B  while  de¬ 
creasing  in  C  over  the  12  hour  period.  The  effect  of  ver¬ 
tical  motion  within  the  volume  is  to  increase  available  poten¬ 
tial  energy  by  increasing  the  strength  of  frontal  zones. 

This  is  seen  in  volume  A,  Figure  5  where  middle  and  upper 
fronts  intensify  while  the  front  at  low  levels  weakens.  In 
volume  C  the  opposite  is  observed  as  the  lower  front  increases 
with  respect  to  a  weakening  upper  front. 
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=£>  300  MB  JET 

v - V~  50  0  MB  FRONT 

T - SFC  FRONT 


JET  MAXIMA 
0)  AT  5000  M 
BUDGET  VOLUMES 


Fig.  3.  Composite  chart  illustrating  surface  and  500  mb  frontal 
position;  budget  Volumes  A,  B,  C;  Tropospheric  jet 
axis  and  wind  maxima;  and  vertical  motion  (mb”3/Sec) 
for  OOZ  6  Nov.  73  (see  key). 


Fig.  4.  Composite  chart  illustrating  surface  and  500  mb  frontal 
position;  budget  volumes  A,  B,  C;  Tropospheric  jet  axis 
and  wind  maxima;  and  vertical  motion  (mb~3/sec)  for  12Z 
6  Nov.  73(see  key). 
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MEAN  QUANTITY 

ooz 

12 Z 

12  HR.  CHANGE 

0)  (MBAR/SEC  x  1(T3) 

-1.4 

-2.7 

-1.3 

KE  (CAL/GRM) 

.088 

.096 

+  .008 

UPPER  HALF 

.123 

.160 

+  .037 

LOWER  HALF 

.053 

.031 

-  .022 

INTERNAL  (CAL/GRM) 

59.740 

58 . 500 

-1.24 

POTENTIAL  (CAL/GRM) 

13.800 

15.000 

+  1.20 

LATENT  HEAT  (CAL/GRM) 

1.13 

1.19 

+  .06 

ABS.  VORT.  (SEC_1  x  ICT5) 

12.5 

11.2 

-1.3 

REL.  VORT.  ( SEC~ 1  x  10~5) 

1.2 

.5 

-  .7 

UPPER  HALF 

-1.3 

2.3 

+3.6 

LOWER  HALF 

3.3 

-1.4 

-4.7 

FRONTAL  INTENSITY  [-|—] 
(°C/111KM) 

3.1 

4  .1 

+  1.0 

UPPER  HALF 

4.0 

5.7 

+1.7 

LOWER  HALF 

2.2 

2.5 

+  .3 

Fig.  5.  Schematic  of  Volume  A  (containing  frontal  zone)  passing 
over  Alps.  Table  lists  mass  weighted  volume  averages 
for  00  and  12Z  6  Nov.  73  and  time  changes  in:  vertical 
motion;  kinetic  energy  (in  upper  and  lower  half  of  volume), 
other  energy  quantities;  absolute  and  relative  vorticity 
(in  upper  and  lower  half);  and  frontal  intensity  (in  upper 
and  lower  half ) . 
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MEAN  QUANTITY 

OOZ 

12Z 

12  HR.  CHANGE 

w  (MBAR/SEC  x  10~3) 

-  .6 

-1.7 

-1.1 

KE  (CAL/GRM) 

.053 

.046 

-  .007 

UPPER  HALF 

.083 

.077 

-  .006 

LOWER  HALF 

.024 

.016 

-  .008 

INTERNAL  (CAL/GRM) 

59.95 

60.08 

+  .13 

POTENTIAL  (CAL/GRM) 

13.47 

14.13 

+  .56 

LATENT  HEAT  (CAL/GRM) 

1.19 

1.02 

-  .17 

ABS .  VORT.  (SEC-1  x  10~5) 

13.5 

11.8 

-1.7 

REL.  VORT.  (SEC_1  x  10~5) 

2.4 

1.9 

-  .5 

UPPER  HALF 

3.0 

3.5 

+  .5 

LOWER  HALF 

1.2 

0.2 

-1.0 

FRONTAL  INTENSITY 

(°C/111  KM)  dy 

3.0 

3.0 

0 

UPPER  HALF 

4.0 

5.8 

+1.8 

LOWER  HALF 

1.9 

.2 

-1.7 

Fig,  6.  Schematic  of  Volume  B  (containing  frontal  zone  passing 

over  Pyrenees,  Table  lists  mass  weighted  volume  averages 
for  00  and  12Z  6  Nov.  73  and  time  changes  in:  vertical 
motion;  kinetic  energy  (in  upper  and  lower  half  of  volume), 
other  energy  quantities;  absolute  and  relative  vorticity 
(in  upper  and  lower  half);  and  frontal  intensity  (in 
upper  and  lower  half). 


15 


MEAN  QUANTITY 

QOZ 

12Z 

12  HR.  CHANGE 

0)  {MBAR/SEC  x  10“3) 

-  .6 

-1.7 

-1.1 

KE  (CAL/GRM) 

.053 

.046 

-  .007 

UPPER  HALF 

.083 

.077 

-  .006 

LOWER  HALF 

.024 

.016 

-  .008 

INTERNAL  (CAL/GRM) 

59.95 

60.08 

+  .13 

POTENTIAL  (CAL/GRM) 

13.47 

14.13 

+  .56 

LATENT  HEAT  (CAL/GRM) 

1.19 

1.02 

-  .17 

ABS.  VORT.  (SEC_1  x  10~5) 

13.5 

11.8 

-1.7 

REL.  VORT.  (SEC”1  x  10*5) 

2.4 

1.9 

-  .5 

UPPER  HALF 

2.8 

6.1 

+  3.3 

LOWER  HALF 

3.7 

2.6 

-1.1 

FRONTAL  INTENSITY  [-—] 
(°C/111KM)  dY 

3.1 

4.5 

+  1.4 

UPPER  HALF 

4.0 

5.0 

+1.0 

LOWER  HALF 

1.9 

4.3 

+  2.4 

Fig.  7. 
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Cross  sections  through  each  volume  are  shown  in  Figures 
8  through  13.  Note  how,  in  B  (Figures  10  and  11),  the  upper 
front  strengthens  while  the  lower  front  appears  to  decouple 
and  weaken.  A  similar  trend  is  seen  in  A  (Figures  8  and  9). 

The  net  loss  of  vorticity  and  potential  vorticity  in 
A  and  B  at  low  levels  appears  to  be  related  to  the  same 
mechanism  (frictional  effect),  which  is  reducing  the  kinetic 
energy.  Evaluation  of  vorticity  tendencies  for  each  data 
time  could  not  isolate  the  term  responsible  for  this  loss,  so 
frictional  dissipation  may  be  the  mechanism. 

Lagrangian  tendencies  of  KE  and  heat  (Figures  14  and 
15)  show  how  the  terrain  influences  the  air  parcels.  Values 
shown  are  for  columns  one  meter  square  and  as  deep  as  the 
analyzed  atmosphere.  Note  how  at  00  and  12Z  the  mountain 
ranges  act  as  KE  sinks.  This  supports  the  results  of  the 
bulk  calculations  in  Figures  5,  6  and  7.  In  addition  the 
forcing  of  air  over  the  mountain  cools  air  on  the  windward 
side  and  warms  air  on  the  leeward  side.  Mountain  ranges  are 
frontogenetic  in  this  respect.  Figures  10  and  11  illustrate 
how  a  weak  low  level  front  is  enhanced  within  B,  as  northerly 
flow  passes  over  the  mountain  barrier. 

Case  2,  21/22  July  1976 

Figures  16  and  17  show  the  synoptic  features  on  21 
July/12Z  and  22  July/OOZ.  The  front  interacts  with  the 
mountain  ranges  in  a  very  similar  way  to  Case  1.  Volumes  A, 

B  and  C  are  assigned  to  appropriate  frontal  elements. 
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10km  — 


Fi».  8.  Cross  section  through  Volume  A  (along  HE)  from  36°  N 

to  55°  N,  illustrating  frontal  zone,  potential  temperature 
(0),  vertical  motion  (w)  and  westerly  wind  component  (u) 
for  00Z  6  Nov.  73. 
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Fig.  9.  Cross  section  through  Volume  A  (along  11  E)  from  36°  N 
to  55°  N,  illustrating  frontal  zone,  potential  temper¬ 
ature  (9),  vertical  motion  (w)  and  westerly  wind  component 
(u)  for  12Z  6  Nov.  73. 
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10km 


Fig.  10.  Cross  section  through  volume  B  (along  1°  W)  from  36°  N 

to  55°  N  illustrating  frontal  zone,  potential  temperature 
(0),  vertical  motion  (go)  and  westerly  wind  component  (u) 
for  00Z  6  Nov .  73. 
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Volume  B  (along  1°  W)  from  36°  N 
frontal  zone,  potential  temperature 
(w)  and  westerly  wind  component  (u) 

or  1 2 Z  6  Nov.  73. 


Fig.  11.  Cross  section  through 
to  55°  N  illustrating 
(9),  vertical  motion 


21 


3  km-  300 


295 

I  km— 


10  km  - 


20m/fcec  30rn/$ec  40nh/feec v 


7km  — 


5  km  — 


30$ 


4  km— 


FRONTAL  SURFACf 
uJ  (  ;j  BAR/SEC) 


30rn/Sec 


i‘ 

tERR^rk! 


Fig .  12. 


Cross  section  through  Volume  C  (along  5°  E) 
55°  N  illustrating  frontal  zone,  potential 
(0)  ,  vertical  motion  (to)  and  westerly  wind 
(u)  for  00Z  6  Nov.  73. 
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<^tl0n  through  Volume  c  {along  5°  E)  from  36°  to 
55  N  illustrating  frontal  zone,  potential  temperature 

f  >,T?Jr*1Sal  I!1^ion  and  westerly  wind  component  (u) 

tor  LIZ  b  Nov.  73  . 
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Fig.  14. 


Lagrangian  kinetic  energy  tendency  through  whole  atmos 
phene  depth.  For  6  Nov.  73.  s 

d (4v2p  )  ...  2  4V 

r  ~  s  (cal/m  sec  x  10  ) 

1  dt  J 
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Fig .  15. 


Lagrangian  enthalpy  tendency  through  whole 
depth  for  6  Nov.  73. 
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Fig.  16.  Composite  chart  illustrating  surface  and  500  mb  frontal 
position;  budget  Volumes  A,  B,  C;  Tropospheric  jet 
axis  and  wind  maxima;  and  vertical  motion  (mb“'^/sec) 
for  12Z  21  July  76. 


Fig.  17.  Composite  chart  illustrating  surface  and  500  mb  frontal 
position;  budget  Volumes  A,  B,  C;  Tropospheric  jet  axis 
and  wind  maxima;  and  vertical  motion  (mb**Vsec)  for  00Z 
22  July  76. 
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Some  similarities  are  observed  (Figures  18  through 
26) .  The  lower  level  kinetic  energy  and  vorticity  is  dimin¬ 
ished  most  significantly  in  mountain  regions.  Frontal  inten- 
®i£y  increases  in  that  part  of  the  front  directly  over  the 
mountain  crest. 

The  dynamical  forcing  induced  by  a  wind  speed  maxima 
in  the  upper  flow  is  clearly  seen  in  Figure  27  as  a  couplet 
of  increasing/decreasing  KE  tendencies  as  air  accelerates 
upstream  of  the  maximum  and  decelerates  downstream.  This  is 
consistent  with  models  of  wind  maxima  (Palmein  and  Newton, 
1969).  This  couplet  moves  across  Europe,  contrasting  with 
the  non-moving  sinks  associated  with  the  mountains.  Figure 
28  shows  the  influence  of  the  mountains  on  heating  and  cool¬ 
ing  of  parcels . 


Conclusion 

While  the  study  completed  is  not  thorough  enough  to 
support  broad  generalizations,  some  useful  findings  did 
result . 

a)  The  analysis  scheme,  in  terms  of  positioning  of  meteoro¬ 
logical  features  (frontal  positions,  jet  positions,  ver¬ 
tical  motion  fields)  seem  to  be  excellent.  The  0)  patterns 
were  quite  good  in  comparison  with  obvious  upper  level 
f rontogenet ic  processes  and  terrain  features.  The  cross 
section  analyses  confirm  this  quite  well.  This  study 
has  increased  confidence  in  the  use  of  basic  analysis 
scheme . 
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MEAN  QUANTITY 

21/12Z 

22/00Z 

12  HR.  CHANGE 

w  (MBAR/SEC  x  10~3) 

-1.5 

-  .9 

+  .6 

KINETIC  (CAL/GRM) 

.023 

.023 

0 

UPPER  HALF 

.030 

.039 

+  .009 

LOWER  HALF 

.015 

.007 

-  .008 

INTERNAL  (CAL/GRM) 

62.66 

61 .56 

-  .90 

POTENTIAL  (CAL/GRM) 

14 .13 

15.33 

1.20 

LATENT  HEAT  (CAL/GRM) 

2.00. 

1.96 

-  .04 

ABS.  VORT. 

10.4 

9.8 

-  .6 

P.EL  VORT.  (SEC’1  x  1CT5) 

-  .7 

-  .6 

+  .1 

UPPER  HALF 

-3.0 

-2 . 3 

-  .7 

LOWER  HALF 

1.5 

1.1 

-  .4 

FRONTAL  INTENSITY  [—-) 

( °C/111KM)  °Y 

2.4 

.9 

-1.5 

UPPER  HALF 

1.3 

.3 

-1.0 

LOWER  HALF 

3.5 

1.5 

-2.0 

Fig.  18.  Schematic  of  Volume  A  (containing  frontal  zone)  passing 
over  Alps.  Table  lists  mass  weighted  volume  averages 
for  1 2 Z  21  July  76  and  00Z  22  July  76  and  time  changes 
in.  vertical  mo  t ion }  kinetic  energy  (in  upper  and  lower 
half  of  volume),  other  energy  quantities;  absolute  and 
relative  vorticity  (in  upper  and  lower  half);  and  fron¬ 
tal  intensity  (in  upper  and  lower  half). 
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MEAN  QUANTITY 

21/12Z 

2  2/00  Z 

12  HR.  CHANG: 

uj  (MBAR/SEC  x  10  ) 

+1.50 

+  .54 

-  .96 

KINETIC  (CAL/GRM) 

.019 

.022 

+  .003 

UPPER  HALF 

.022 

.033 

+  .011 

LOWER  HALF 

.015 

.010 

-  .005 

INTERNAL  (CAL/GRM) 

62.50 

61.49 

-1.01 

POTENTIAL  (CAL/GRM) 

14.50 

14 .74 

+  .24 

LATENT  HEAT  (CAL/GRM) 

1.80 

1.72 

-  .08 

ABS.  VORT.  (x  10~5SEC_1) 

12.9 

9.4 

3.5 

REL.  VORT.  (x  10-5SEC_1) 

+1.5 

-  .5 

-2.0 

UPPER  HALF 

+  2.0 

0 

-2.0 

LOWER  HALF 

+1.0 

-  .9 

-1.9 

FRONTAL  INTENSITY 

{ °C/111KM)  dy 

2.2 

1.5 

-  .7 

UPPER  HALF 

.5 

1.0 

+  .5 

LOWER  HALF 

1.5 

2.0 

+  .5 

Schematic  of  Volume  B  (containing  frontal  zone  passing 
over  Pyrenees) .  Table  lists  mass  weighted  volume 
averages  for  12Z  21  July  76  and  00Z  22  July  76  and  time 
changes  in:  vertical  motion;  kinetic  energy  (in  upper 
and  lower  half  of  volume),  other  energy  quantities; 
absolute  and  relative  vorticity  (in  upper  and  lower  half 
and  frontal  intensity  (in  upper  and  lower  half). 
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MEAN  QUANTITY 

21/22  Z 

2  2/00  Z 

12  HR.  CHANGE 

o)  (MBAR/SEC  x  10~3) 

+  .2 

+1.7 

+1.5 

KINETIC  (CAL/GRM) 

.028 

.024 

-  .004 

UPPER  HALF 

.043 

.039 

-  .004 

LOWER  HALF 

.011 

.010 

-  .001 

INTERNAL  { CAL/GRM ) 

61.88 

61.78 

-  .10 

POTENTIAL  (CAL/GRM) 

14 .72 

14 .71 

-  .01 

LATENT  HEAT  (CAL/GRM) 

1.41 

1.51 

+1.0 

ABS.  VORT.  (SEC'1  x  10~5) 

13.7 

12.0 

-1.7 

REL.  VORT.  (SEC"1  x  10”5) 

3.0 

3.0 

0 

UPPER  HALF 

4.5 

4.0 

-0.5 

LOWER  HALF 

1.5 

2.0 

+  .5 

FRONTAL  INTENSITY 

( °C/1 11KM)  dy 

2.0 

3.2 

+1.2 

UPPER  HALF 

2.0 

4 .0 

+2.0 

LOWER  HALF 

2.0 

2.5 

+  .5 

Fig.  20.  Schematic  of  Volume  C  moving  into  Mediterranean.  Table 
lists  mass  weighted  volume  averages  for  12Z  21  July  76 
and  00Z  22  July  76  and  time  changes  in:  vertical  motion; 
kinetic  energy  (in  upper  and  lower  half  of  volume),  other 
energy  quantities;  absolute  and  relative  vorticity  (in 
upper  and  lower  half);  and  frontal  intensity  (in  upper 
and  lower  half). 
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U)  (  f, bar/sec  ) 


Fig.  21.  Cross  section  through  Volume  A  (along  11E)  from  36°  N 
to  55°  N,  illustrating  frontal  zone,  potential  temper¬ 
ature  (0),  vertical  motion  (w)  and  westerly  wind  com¬ 
ponent  (u)  for  12Z  21  July  76. 
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10  km 


4  km— 


Fig.  22. 


Cross  section  through  Volume  A 
to  55°  N,  illustrating  frontal 
ature  (6),  vertical  motion  (u>) 
(u)  for  00Z  22  July  76. 


(along  HE)  from  36°  N 
zone,  potential  temper- 
and  westerly  wind  component 
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JO  km  — 


Fig.  23.. 


Cross  section  through 
to  55"  K  illustrating 
(0)  ,  vertical  motion 
for  12Z  21  July  76. 


Voiutne  B  (along  1°  w)  from  36"  N 
frontal  z one,  potential  tempers  tut 
(u)  and  westerly  wind  component  (u) 
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Fig,  24.  Cross  section  through  Volume  B  (along  1°  W)  from  36°  N 

to  55°  N  illustrating  frontal  zone,  potential  temperature 
(0),  vertical  motion  (go)  and  westerly  wind  component  (u) 
for  00Z  22  July  76. 
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Fig.  25.  Cross  section  through  Volume  C  (along  5°  E)  from  36°  N 

to  55  N  illustrating  frontal  zone,  potential  temperature 
(0),  vertical  motion  (u)  and  westerly  wind  component  (u) 
for  1 2 Z  21  July  76. 
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Fig.  26.  Cross  section  through  Volume  C  (along  5°  E)  from  36°  N 
to  55°  N  illustrating  frontal  zone,  potential  tempera¬ 
ture  (6),  vertical  motion  (to)  and  westerly  wind  com¬ 
ponent  (u)  for  00Z  22  July  76. 
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Fig'.  27.  Lagrangian  kinetic  energy  tendency  through  whole  atmos¬ 
pheric  depth.  For  21-22  July  76. 

d(yv2P  )  9  , 

[ - - ]  (cal/m  sec  x  10  ; 
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Fig.  28.  Lagrangian  enthalpy  tendency  through  whole  atmospheric 
depth  for  21-22  July  76. 

d(C  TP  ) 

[  ^ — ~]  (cal/m  sec  x  10  ) 


38 


b)  Energy  budgets,  while  not  balancing  perfectly,  contained 
only  small  errors  which  in  most  cases  could  be  explained 

by  energy  source  terms  such  as  daytime  heating  or  friction. 

c)  Kinetic  energy  reduction  in  the  lower  half  of  each  volume 
compares  closely  with  the  "roughness"  or  maximum  height 
of  the  terrain  encountered.  See  Figure  29. 

d)  Not  surprisingly,  calculations  of  parcel  tendencies  in¬ 
dicate  mountains  act  as  KE  sinks,  and  also  as  heat  sinks 
on  the  windward  side  and  heat  sources  in  the  lee.  The 
latter  phenomena  tends  to  enhance  the  baroclinic  zone  or 
front  as  it  passes  over  the  mountain  ridge.  Figure  30 
shows  a  schematic  diagram  of  this  mechanism. 

e)  The  net  effect  of  moisture  condensation  in  a  frontal  zone 
passing  over  terrain  (higher  than  the  lifted  condensation 
level)  is  to  slow,  steepen  and  intensify  the  front. 

Figure  31  shows  in  a  simple  way  how  this  is  done.  Com¬ 
parison  of  the  November  (dry)  case  and  the  July  (moist) 
case  points  this  out.  Fronts  weakened  less  (even  inten¬ 
sified)  in  low  levels  after  the  front  passed  the  mountain 
barrier  in  the  July  case. 

The  moisture  effect  may  also  be  seen  in  Case  2, 
volume  A  where  the  windshift  line  moved  far  in  advance 
of  the  temperature  gradient.  The  destabilization  of  air 
ahead  of  the  front  is  a  commonly  seen  feature  in  Pacific 
fronts,  often  producing  strong  gusty  winds  east  of  the 
Rockies  in  the  dry  air  preceding  the  front. 
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TERRAIN  HEIGHT  (  M  ) 


Fig.  29.  Percent  reduction  (for  12  hours)  in  bulk  kinetic 

energy  in  volumes  A,  B,  C  for  Case  1  (Nov.  73)  and 
Case  2  (Jul.  76)  versus  maximum  terrain  height 
along  volume  trajectory.  (Plot  A  indicates  Volume 
A  for  Case  1)  . 
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Schematic  diagram  illustrating  the  effect 
tain  ridge  on  frontal  intensity.  Such  an 
seen  in  the  case  study  analysis. 


of  a  moun- 
effect  was 
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FRONTAL  POSITION 

Fig.  31.  Schematic  diagram  illustrating  the  effect  of  moisture 
on  frontal  intensity  r  3  6-.  passing  a  mountain  ridge. 
The  effect  was  seen  3yJ  in  the  case  study  analysis. 
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f)  The  use  of  vorticity  equations  as  a  diagnostic  tool  was 
not  too  successful.  Increases  in  vorticity  above  the 
mountain  crest  observed  in  some  cases  were  partially  ex¬ 
plained  by  the  tilting  term,  but  continuity  with  respect 
to  potential  vorticity  was  difficult  to  interpret.  Loss 

of  low  level  vorticity  may  be  due  to  mountain  frictional 
effects • 

g)  Heat  and  kinetic  energy  flux  data  were  gathered,  but  time 
limits  prevented  meaningful  interpretation. 
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APPENDIX  A 


USE  OF  VARIATIONAL  INEQUALITY  CONSTRAINTS  IN 
ADJUSTMENT  OF  VERTICAL  TEMPERATURE  PROFILES 

The  goal  of  this  work  in  objective  analysis  is  to 
impose  upper  o_r  lower  limits  on  one  or  more  meteorological 
quantities,  this  being  combined  with  observational  and  dy¬ 
namical  constraints  which  have  been  well  documented  (Sasaki, 
1958  ,  1969  ,  1970)  . 

This  class  of  constraint  has  been  known  and  in  use 
in  engineering,  particularly  in  control  problems. 

The  strong  inequality  constraint  can  be  reduced  to  a 
more  familiar  strong  constraint  problem  or  Bolza  problem  by 
the  use  of  the  Val en t ine-Berkovit z  Theorem,  which  introduces 
an  unknown  function  termed  the  "slack"  function  (Berkovitz, 
1961) . 

To  illustrate,  consider  the  present  problem.  Temper¬ 
ature  fields  are  objectively  analyzed  using  low  pass  varia¬ 
tional  filters  (Wagner,  1971).  However,  in  the  vertical, 
instrument  problems  or  low  data  density  occasionally  produce 
spurious  super  adiabatic  layers.  For  a  variety  of  reasons 
it  is  desirable  to  remove  these  layers  by  re-adjusting  the 
temperatures  to  conform  to  neutral  stability  constraints 
(which  are  consistent  with  the  scale  of  our  analysis). 
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Thus  the  condition  on  the  temperature  becomes: 


8  it  ^  ®  where  0  is  potential  temperature  and  ir  is  the 
vertical  coordinate. 

At  this  point  we  introduce  the  unknown  slack  function 
D(tt),  which  allows  us  to  write  the  constraint  as  an  equality: 


g  0 

§T  =  d2^)  N°te  that  D 2  (tt)  (if  real  valued)  is  posi- 
tive  for  all  tt  . 


The  variational  formalism  can  now  be  expressed  as: 


TT 


w 

J7 


{ ot  ( tt )  (0-0)  2  +  X  ( 7r )  (|^  -  D2  (ir)  )  ]  dir 


8  TT 


where  a  is  an  observational  weight  function  (a  function  of 
U)  »  6  ls  the  "observed"  or  unadjusted  temperature,  X(tt)  is 
the  Lagrange  multiplier  associated  with  the  strong  constraint 

stated  above.  0,  X,  and  D  and  all  unknown  functions  which 
must  be  determined. 


Taking  the  variation  and  setting  it  equal  to  0  allows 
us  to  find  the  functions  X,  D,  0  which  minimize  J. 


TT 


6  J  = 


/ 


[2  (0-0)  60  +  (||  -  D2)  6X  +  X  ^  -  2 X D 6 D ]dir  =  0 


After  integrating  by  parts  and  collecting  terms,  and  assuming 
the  variations  are  arbitrary  through  the  domain  [  0 , 7TQ  ]  #  (we 
assume  natural  boundary  conditions  on  X  [  \  =  X  =  o])  the 
three  Eul er-Lagrange  equations  are  obtained. 
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1)  2a  (9  -  0) 


0 


3_\ 
3  TT 


2) 


||  -  D2 

3  IT 


=  0 


3)  Ad  =  0 


Combining,  we  get: 


3  2X 

3  TT2 


-  2a (D2  -  |f)  - 


2  |f  (9  -  8)  -  0 


This  equation  is  put  into  finite  difference  form  and  solved 
by  an  iterative-relaxation.  Throughout  the  computation  the 
condition  expressed  by  eqn.  3  is  imposed. 

Figure  A1  shows  unadjusted  and  adjusted  0  profiles. 

The  weight  a  can  be  set  for  each  level  depending  on  the 
quality  of  data.  Figure  A3  shows  with  a  low  weight 

(suspected  instrument  error:  possible  wet  temperature 

element  at  7^  )  .  Figure  A2  shows  0(7^)  with  a  weight  equal 
to  other  levels  for  circumstances  where  the  atmosphere  (or 
analysis  scheme)  may  be  producing  short  lived  super  adiabatic 
layers . 

It  should  also  be  noted  that  if  a(7T)  is  constant  for 

all  7T  then  the  condition  A.  =  X  =  0  forces  conservation  of 

0  Tf  o 

the  mean  potential  temperature  within  the  layer. 
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Fig.  Al.  Hypothetical  "observed"  and  adjusted  potential  tem¬ 
perature  profiles  using  inequality  constraint. 
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<X(6)  = 


7l 


oc(5)  = 
*(4)  = 
«(3)  = 
<X(2)  = 
0C(  I  )  = 
ot(  0 )  = 


. -  OBSERVED  9  PROFILE 

-  ADJUSTED  6  PROFILE 

Fig.  A2.  Hypothetical  "observed"  and  adjusted  potential  tem¬ 
perature  profiles  for  equal  weighting  case  (a(K)  for 
K  =  1,6).  Here  data  at  3  is  assumed  to  be  due  to 
actual  (though  possibly  short  lived)  atmospheric  con¬ 
ditions  . 


0t(4)  =  I 
0C(3)  =  .1 
0t(  2)  =  I 


Fig.  A3.  Hypothetical  "observed"  and  adjusted  potential  tempera¬ 
ture  profiles  for  non  equal  weighting  ( ct  ( 3 )  =  .1 
c*(K)  =  1.0)  where  error  may  be  due  to  systematic  instru¬ 
ment  error  such  as  rawinsonde  exiting  cloud  layer. 
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